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ABSTRACT 

Association of the high-frequency quasi-periodic oscillation (HFQPO) pairs with the 
broad Fe K line in XTE J1550-564 and GRO J1655-40 is discussed based on the 
magnetic coupling (MC) of a rotating black hole (BH) with its surrounding disc. The 
3:2 HFQPO pairs are interpreted by virtue of the inner and outer hotspots arising from 
non-axisymmetric magnetic field, where the inner hotspot is produced by a torque 
exerted at the inner edge of the disc, and the outer hotspot is created by the screw 
instability of the large-scale magnetic field. The very steep emissivity index is created 
predominantly by the torque exerted at the inner edge of the disc. It turns out that 
the 3:2 HFQPO pairs observed in the two sources can be fitted by tuning several 
model parameters, such as the BH spin, and the main features of this model lie in 
three aspects. (1) The condition for only one HFQPO is discussed based on the two 
mechanisms for producing the 3:2 HFQPO pairs, (2) an explanation is given for a 
systematic shift away from disc dominated flux with the increasing power-law flux as 
the HFQPO pairs shift from the higher to lower frequencies, which is consistent with 
the analysis given by Remillard et al. (2002), and (3) the BH spin in XTE J1550-564 
and GRO J1655-40 can be estimated by combining the 3:2 HFQPO pairs with the 
very steep emissivity index required for fitting the broad Fe K emission line. 

Key words: accretion, accretion discs — black hole physics — magnetic fields — 
stars: individual (XTE J1550-564, GRO J1655-40) — - stars: oscillations ~ X-rays: 
stars 



1 INTRODUCTION 

Quasi-periodic oscillations in black hole X-ray binaries (BHXBs) have become a very active research field since the launch of 
the NASA satellite RXTE (Bradt et al. 1993). One of the remarkable features in BHXBs is that the high-frequency quasi- 
periodic oscillations (HFQPOs) could appear in pairs with the puzzhng commensurate frequencies, which have been observed 
in GRO J1655-40 (450, 300Hz; Remillard et al. 1999; Strohmayer 2001a; Remillard et al. 2002, hereafter R02), XTE J1550-564 
(276, 184, 92Hz; Miller et al. 2001; R02) and GRS 1915+105 (168, 113, 67, 41Hz; McClintock & Remillard 2006; Remillard 
& McClintock 2006, hereafter MR06 and RM06, respectively). The above HFQPOs appear in pairs with the 3:2 ratio of 
the higher frequency to the lower frequency (henceforth 3:2 HFQPO pairs), e.g. 276 vs 184Hz and 450 vs 300Hz occur in 
XTE J1550-564 and GRO J1655-40, respectively, while two 3:2 HFQPO pairs, 168 vs 113Hz and 67 vs 41Hz, appear in GRS 
1915-hl05. 

A number of models have been proposed to explain the origin of HFQPO pairs in BHXBs. Strohmayer (2001a, 2001b) 
investigated combinations of the azimuthal and radial coordinate frequencies in general relativity to explain the HFQPO pairs 
in GRO J1655-40 and GRS 19154-105. Wagoner et al. (2001) regarded the HFQPO pairs as fundamental g-mode and c-mode 
discoseismic oscillations in a relativistic accretion disc. Very recently, Silbergleit & Wagoner (2007) discussed the corotation 
resonance and diskoseismology modes of black hole (BH) accretion discs and suggested that the HFQPO pairs may relate 
to the excitation of two (groups of) g-modes of discoseismic oscillations. Abramowicz & Kluzniak (2001) explained the pairs 
in GRO J1655-40 as a resonance between orbital and epicyclic motion of accreting matter. Recently, the resonance model is 
presented in a more realistic context, in which "parametric resonance" concept is introduced to describe the oscillations rooted 
in fluid flow where there is a coupling between the radial and polar coordinate frequencies (Abramowicz et al. 2003; Kluzniak 
et al. 2004; Torok et al. 2005). As argued by van der Klis (2000, 2006), HFQPOs in X-ray binaries probably originate from 
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the inner edge of an accretion disc around a neutron star or a stellar-mass BH, since millisecond is the natural timescale for 
accretion process in these regions. 

On the other hand, Wilms et al. (2001) presented the first XMM-Newton observation of MCG-6-30-15, and they found a 
broad Fe disc line from the observed spectrum. However, no quasi-periodic oscillations have been detected in this source. A 
very stoop omissivity index (4.3 < /3 < 5.0) from the irmcr accretion disc is required for fitting the broad Fc K omission line, 
which is difficult to be explained within the framework of a standard accretion disc (SAD). Wilms et al. (2001) suggested that 
the magnetic extraction of rotational energy from a spinning black hole should be invoked to create this steep emissivity. 

The relativistic disc lines for microquasaxs have been analyzed by some authors, e.g., XTE J1550-564 by Sobczak et al. 
(2000), Miller et al. (2003, 2004) ; XTE J1650-500 by Miller et al. (2004), Miniutti, Fabian, and Miller (2004); GRO J1655-40 
by Miller et al. (2004), Diaz Trigo ct al. (2007). However, the steep emissivity indexes required by broad Fc K a lines were 
hardly worked out except a few cases, such as, (3 = 5.5^0.7 for GX 339-4 in the very high state, and /3 ~ 5 for XTE J1650-500 
based on the observation of XMM-Newton. A detailed review for relativistic X-ray emission lines from the inner accretion 
disc around black holes for AGNs and stcllar-mass black holes arc given by Miller (2007). 

With the existence of a magnetic field connecting a rotating black hole to its surrounding disc, energy and angular 
momentum can be transferred from the black hole to the disc, and this energy mechanism is referred to as the MC process, 
which is regarded as one of the variants of the Blandford-Znajek (BZ) process (Blandford & Znajek 1977, 1999; Li 2000, 2002, 
hereafter L02, Wang 2002). 

Wang et al. (2002, 2003a, hereafter W03a) incorporated the BZ and MC processes into black hole accretion disc, and 
henceforth this model is referred to as MC-I model, which was used to interpret the very steep emissivity index required by the 
observation of MCG-6-30-15. In addition, Wang et al. (2005, hereafter W05) suggested that the 3:2 HFQPO pairs observed 
in the above BHXBs can be fitted by virtue of two hotspots in the inner disc, which are produced by the non-axisymmetric 
MC process. 

The fatal shortcoming of MC-I model is that the fits of the 3:2 HFQPO pairs always accompany jots, which is driven by 
the BZ process. As reviewed in MR06 and RM06, jets are generally observed in the hard state of BHXBs, while HFQPOs 
are always detected in the steep power-law (SPL) state. Thus the fits of the 3:2 HFQPO pairs based on MC-I model are 
inconsistent with the observations. 

Based on a careful analysis of the 1998-1999 outburst of XTE J1550-564 and the 1996-1997 outburst of GRO J1655-40 
some interesting features have been found in R02, the 3:2 HFQPO pairs could bo dotoctod simultaneously in the SPL state, 
and a systematic shift away from disc dominated flux with the increasing power-law fiux is detected as the HFQPO pairs 
shift from 276 to 184 Hz for XTE J1550-564 and from 450 to 300 Hz for GRO J1655-40 as shown in Figs. 5 and 7 of R02, 
respectively. In addition, Remillaxd (2004) noted that the broad Fe K emission lines can be also detected in the SPL state in 
XTE J1550-564 and GRO J1655-40. These results imply that the Fe K lines could be detected simultaneously with the 3:2 
HFQPO pairs in these two sources. 

In this paper, we interpret the 3:2 HFQPO pairs associated with the broad Fe K line in XTE J1550-564 and GRO 

J1655-40 by modifying MC-I model in two aspects. (1) The open field lines corresponding to the BZ process is removed, 
and (2) a non-zero torque exerted at the inner edge of the disc is introduced. Henceforth the modified model is referred to 
as MC-II model. It turns out that the 3:2 HFQPO pairs observed in the two sources can be fitted by tuning several model 
parameters, such as the BH spin, and the main features of MC-II model lie in the following aspects: (1) The condition for 
only one HFQPO is discussed based on the two mechanisms for producing the 3:2 HFQPO pairs. (2) An explanation is given 
for a systematic shift away from disc dominated fiux with the increasing power-law flux as the HFQPO pairs shift from the 
higher to lower frequencies, which is consistent with the analysis given by R02. (3) The BH spin in XTE J 1550-564 and GRO 
J1655-40 can be estimated by combining the 3:2 HFQPO pairs with the very steep emissivity index required for fitting the 
broad Fe K emission line. 

This paper is organized as follows. In §2 we present a detailed description of MC-II model, in which the main differences 
between MC-I and MC-II models are outlined. In §3 we describe the elements for fitting the 3:2 HFQPO pairs observed in 
XTE J1550-564 and GRO J1655-40, and present the calculation sequence and results based on MC-II model. We discuss the 
possible correlation of the 3:2 HFQPO pairs with the very steep emissivity index required by broad Fe K a lines. It is expected 
that the 3:2 HFQPO pairs could be observed with the steep emissivity index, and the black hole spin could be estimated 
once the value ranges of the index are determined in the future observations. Finally, we discuss some issues related to MC-II 
model in §4. 

Throughout this paper the geometric units G = c = 1 are used. 



2 DESCRIPTION OF MC-II MODEL 

In this section we present a detailed description of MC-II model based on the following assumptions. 

(1) The accretion disc is perfectly conducting, and the closed magnetic field lines are frozen in the disc. The disc is thin 
and relativistic, lying in the equatorial plane of a Kerr black hole. 

(2) Corona is introduced, which consists of tenuous hot plasma above the disc surface as shown in Figure 1. Both large- and 
small-scale magnetic fields are involved in MC-II model, in which the large-scale magnetic field plays a key role in transferring 
energy and angular momentum from a fast-rotating black hole to its surrounding disc in the MC process, while the small-scale 
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magnetic field energizes the corona above the disc by virtue of some physical processes, such as magnetic reconnection and 
buoyancy (Haardt & Maraschi 1991). It has been shown that the existence of corona can improve the fitting of the output 
spectra from the BHXBs, which is helpful to interpret the X-ray radiation observed in the SPL state (Ma et al. 2006). 

(3) The large-scale magnetic field is assumed to be non-axisymmetric as described by Wang et al. (2003b, hereafter 
W03b). In addition, the magnetic field is assumed to be constant on the black hole horizon, and it varies as a power law with 
disc radius, i.e., Bd oc r~". 

(4) One of the main differences between MC-I and MC-II models lies in the fact that "no torque boundary condition" is 
assumed in MC-I model as given in SAD (Novikov & Thorne 1973, hereafter NT73; Page & Thorne 1974), while a nonzero 
torque is exerted at the inner edge of the disc in MC-II model. This torque arises from the magnetic connection of the matter 
in the plunging region with the matter in the inner disc as argued by some authors (Livio, Ogilivie & Pringle 1999, hereafter 
L99; Krolik 1999; Gammie 1999; Agol & Krolik 2000). 

(5) Another difference between the two models lies in the fact that the open field lines are removed from MC-II model, 
while the closed field lines connecting the black hole with its surrounding accretion disc are retained, as shown in Figure 1. 
So only the MC process rather than the BZ process can work in MC-II model. 

The reason for the modification of the magnetic field configurations lies in two points. 

(i) The existence of the closed field lines depends upon the mapping relation with the conservation of the magnetic fiux 
as argued in W03a, being not the sufficient condition for the existence of the open field lines at the horizon. 

(ii) The second reason is more important, i.e., the 3:2 HFQPO pairs are observed in the SPL state, in which no jets are 
observed. Thus the magnetic field configuration depicted in Figure 1 is more consistent with the observations. 

In order to illustrate the elements for interpreting the 3:2 HFQPO pairs we derive the radiation flux from the accretion 
disc based on MC-II model. According to the conservation of energy and angular momentum we have the following relations 
given in L02 and W03a, 

^iMoL'' - 3„,,) = 47rr(Ft,t,,Lt - Hmc), 

(1) 

^{MdE'' - g^^s^o) = A-Kr{FtotaiE'^ - Hmc^d). 

In equation ((l| Mo is the accretion rate. The quantities g^is and Hmc are respectively the interior viscous torque in the 
disc and the fiux of angular momentum transferred between the black hole and the disc in the MC process, and and 
are the specific energy and angular momentum of the disc matter, being given in NT73 as follows (where x = ^/r/M), 

E' = il- 2x'' + a,X~') /(I - 3x-' + 2a..x~'Y'\ (2) 

= Mx(l - 2a,x-' + «'x"*) /(l - ^x'^ + 2a*x-')'^'. (3) 

By resolving equation ((TJ we derive the total radiation fiux Ftotai and the interior viscous torque g^is in the disc as 
follows, 

Ftotai = Fda + Fmc + Fgi„, (4) 
E^ - »oLt ^ 

gvis = r=; — i-^^Trrl'total- (5) 

In equation Fda and Fmc are respectively the radiation fiuxes contributed by the disc accretion and the MC process, 
and Fgin is the extra contribution due to the exterior torque exerted at the inner edge, and S7_d is the angular velocity of the 
relativistic thin disc and it reads 

= Mix'+a,y 

In equation (|6} the parameter a* = J / is the black hole spin, which is defined in terms of the black hole mass M and 
angular momentum J, and the parameter x is related to the disc radius by x = V^Xms, where ^ = r/rms is defined in terms 
of the innermost stable circular orbit (ISCO, NT73). The three parts of Ftotai are expressed as follows, 

= 1 -dQn/dr r _ ^^^t) tM^dLydr) dr, (7) 



Fmc = -^{E^ -^dL^) ^ I {E^ -nnL^)HMcrdr, (8) 



ir ^ 

1 —dQo/dr / I 



where fli„ = ^Inirin) is the angular velocity at ri„, other terms with subscripts "in" are taken their values at Vin- 

The function Hmc in equation ((8)1 is related to the MC torque by OTmc /'^'^ ~ A-ktHmc, and Tmc is the MC torque due 
to non-axisymmetric magnetic field and it reads (W03b) 
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TmS = Ar,1c, (10) 

where T^iq is the MC torque due to axisymmetric magnetic field, being expressed in W03a as 

r,Vro^4..(i + .) f Os<e<e.. (ii) 

7ss 2- (l-g)sm^6i 

The parameter A in equation (|10|l is used to indicate the difference between T^^^ and T^iq as described in W03b, and 
Pn = ^d/^h is the ratio of the angular velocity of the disc to that of the black hole. It is easy to check by using equation 
([6| that /3n < 1 always holds for a fast-spinning black hole with a, > 0.3594, implying that energy and angular momentum 
are transferred from the black hole to the inner disc in the MC process. 

Since the magnetic field on the black hole is supported by the surrounding disc, there is some relation between Bh and 
Md- One possibility has been suggested by Moderski et al. (1997), which is based upon the balance between the pressure of 
the magnetic field on the horizon and the ram pressure of the innermost parts of an accretion flow, i.e., 

Bl /(Stt) = Pram ~ pc ~ Mo/i'^-nrl) , (12) 



where th is the radius of the black hole horizon. Considering that equation (|12p is not a certain relation between Bh and 
Md, we rewrite it as follows, 

Mn = a^Blrl, (13) 

where am is a parameter to adjust the accretion rate Mn- 

Inspecting equation ([9)l, we find that the contribution Fgin is related directly to the exterior torque gin. According to 
L99 the strength of the magnetic field produced by dynamo process in the disc is given by 



Blyn W 



(14) 



47r 2hin 

and the torque gin is related to the viscous stress at ri„ as follows, 

where hin is the half-thickness at ri„, being related to ri„ by hin = anrin. 

In equation (fH)) W is the integrated shear stress of the disc, and the disc magnetic field -Bij is related to Bdyn by 

Sl3~(V0„.ax5d,„, (16) 

where h is the half-thickness of the disc. As argued in W03a, the field -Bd is related to Bh by 

BD=BH{rH/r^in), (17) 

where LUin is the cylindrical radius at Tin in the context of Kerr metric. 
Incorporating equations (|16|) and p7|) . we have 

Bdyn = {rin/hin)^ Bo = [{rinV h) / {^inhin)]'^ Bjj. (18) 

Incorporating equations (|14|) . (|15|) and psp . we have the torque gin as follows, 
gin = anBdynrfn = ajj^rfn (rn/r^inf B%. (19) 

Assume that the inner edge of the disc is initially located at ISCO. If the magnetic pressure inside ISCO is strong enough, 
a radial force can be exerted at the inner edge of the disc due to the magnetic connection, resulting in an outward displacement 
of the inner edge from Vms to Tin- It is assumed that the magnetic torque gin exerted at Tin is determined by the following 
equations, 

5gin^in = [E'^ (ri„,a*) - E'^ (r„,,a.)] Md, (20) 
where gin^in is the power of the magnetic torque, and S is the fraction of the power keeping the inner edge at ri„. 



3 FITTING 3:2 HFQPO PAIRS ASSOCIATED WITH BROAD EE K LINE 

Considering the fact that the accretion disc is perfectly conducting, and the closed magnetic field lines are frozen in the disc, 
and the inner and outer hotspots are produced by non-axisymmetric magnetic field, we infer that these hotspots rotate with 
the angular velocity of the disc, and the upper and lower frequencies of the 3:2 HFQPO pairs can be expressed by 

vhfqpo = vo {f^^Xms +a*) \ (21) 
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where vq = {mBH)~^ x 3.23 x 10*//^, and tubh = M / Mq is the black hole mass in terms of one solar mass. Equation (|21l) 
can be derived directly from equation (|6]). 

The lower frequency i^iomer of the HFQPO pair is produced by the outer hotspot, which is located at the outer boundary 
of the MC region as shown in Figure 1. According to the Kruskal-Shafranov criterion, the screw instability of the large-scale 
magnetic field will occur, if the magnetic field line turns around itself about once (Kadomtsev 1966; Bateman 1978). Based 
on the Kruskal-Shafranov criterion the radius Vout can be determined by the following equation (Wang et al. 2004), 

(27TmD/L)B^^/BZ ^l, (22) 

where and B]j are the poloidal and toroidal components of the magnetic field on the disc, respectively. The quantities L 
and WD are respectively the poloidal length of the field line and the cylindrical radius on the disc, and the latter reads 

WD = So/pD = ^MxL Vl + aU-'^X^i + 2a2^-3xm., (23) 

where E_d and pD are the Kerr metric parameters (Novikov & Thorns 1973). 

Wang et al. (2004) expressed equation (|22p in terms of the parameters involved in the MC process as follows. 



{2nwD/L)Fsc{a,,^,n) i^l. (24) 



The equality in equation (I24|l corresponds to the screw instability, and Fsc{a*,£,,n) is a function of the parameters, a,, 
^ and n, and it reads 



J, , e-"{l + l) [2csc^g-(l-g)] l + aix^U-^ + 2aix^U-' 



where q = -^/l — is a function of the black hole spin. Thus the position of the outer hotspot can be determined by the 
criterion (|24|) for the given a, and n. 

The upper frequency i^higher of the HFQPO pair is produced by the inner hotspot, of which the position is determined 
by the maximum of the following function, 

FHFQPO=fFtotal/Fo, (26) 

where Fq = Bjj = Bl X ^3 X lO^^erg ■ s ^ ■ cm is defined as a unit of the radiation fiux, and _Bg is the magnetic field in 
terms of lO^gauss. 

The emissivity index can be calculated from the total radiation flux Ftotai based on the deflnition as follows, 

f5^-d\nFtotai/d\nr. (27) 

Incorporating equations Q, ((7|), (O, (O, p9|) and (|27|) . we obtain the emissivity index /3 = Pda + Pmc + Pgin by 
calculating the following contributions, 

_ r dFoA a - AFmc n - ^ dFgi„ 

PDA = -— 1 — , Pmc = -— ^ , I3gin = -— 1 — , (28) 

Ftotai dr Ftotai dr Ftotai dr 

where Pda, Pmc and [3gi„ axe contributed by the disc accretion, the MC process and the magnetic torque exerted at the 
inner edge of the disc, respectively. 

Based on equations (O, ((Si, ([9]) and (|28|) we obtain the curves of the above contributions to the total radiation flux and 
the emissivity index versus the disc radius for XTE J1550-564 with rriBH = 10.8 as shown in Figure 2a and 2b, respectively. 

As shown in Figure 2 the contributions of the disc accretion and the MC process to the total radiation flux and the 
emissivity index are more than two orders of magnitude less than those due to the magnetic torque at the inner edge of the 
disc, and both the total radiation fiux Ftotai and the emissivity index (3 take their maxima at Vin- Thus we infer that both 
the inner hotspot and the maximum emissivity index are located at Tin. 

Combining equation ([21} with the higher and lower HFQPO frequencies corresponding to the disc radii = (.inrms and 

Tout = ioutTms, WC haVC 

yhigher = ^oi^^in xLs + a*)~^ (29) 

yionisr = fo(C!{?xL + a,)~^ (30) 

Equations (|20p . (|24| l. (l27l) . (|29p and H30p can be regarded as a set of equations for calculating the disc radii, ri„ and Tout , 
for the 3:2 HFQPO pairs, and henceforth the five independent equations are referred to as FIE. Four input parameters, a,, 
nfiBH, i^higher aud viomer , are iuvolvcd, aud five output parameters, n, Tin, rout, 5 and /3 are obtained based on FIE. The 
detailed calculation sequence is described as follows. 

(1) As the first step, we assume a value of spin with the observed black hole mass, e.g., a* = 0.5 and m,BH = 10.8 for 
XTE J1550-564. Based on the higher frequency fhigher, e.g., 276Hz for XTE J1550-564, the position of the inner hotspot can 
be determined by equation (|29p . Thus the radius can be determined for the given black hole spin and mass. 

(2) In the second step, the position of the outer hotspot corresponding to the lower frequency uiomer , e.g., 184Hz for XTE 
J1550-564, can be determined by combining equation (|30p with the criterion of the screw instability given by equations (|24p . 
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and the power law index n indicating the variation of the magnetic field with the disc radius can be determined in calculating 
viomer by virtuc of the criterion of the screw instability. 

(3) Finally, in the third step, the outward displacement of the inner edge from rms to ri„ can be worked out, and the 
parameter S is determined by equation (I20|l . Combining equation (|27p with equation and the obtained values a,, Tin and 
n, we have a majcimum value of the emissivity index /3 located at ri„. 

In this way, we can obtain one set of solution of ri„, rout, n, 5 and /3 by inputting each black hole spin a* with the 
observed Vhigher, Slower and tubh into FIE, and the curves of {3, n and 5 varying with a, for XTE J1550-564 and GRO 
J1655-40 are shown in Figures 3 and 4, respectively. 

By resolving FIE we have the specific solutions corresponding to the two BHXBs with the lower and higher black hole 
masses and different spins as listed in Table 1. 

Inspecting Figures 3 and 4, and Table 1, we find that the 3:2 HFQPO pairs observed in XTE J1550-564 and GRO 
J1655-40 can be worked out with the steep emissivity index, and the features of MC-II model are summarized as follows. 

(1) For a given black hole mass, a bigger spin corresponds to a smaller emissivity index (3 and a smaller parameter n, but 
to a bigger parameter 5. For a given black hole spin, a bigger black hole mass corresponds to a bigger 13, while the parameters 
n and 5 are insensitive to the variation of the black hole mass. This result can be understood based on black hole physics 
(Shapiro & Teukolsky 1983). The width of the plunging region decreases monotonically with the increasing black hole spin, 
so the effects of the magnetic torque in the plunging region decreases with the increasing black hole spin. 

(2) The black hole spin could be estimated, if the value range of the emissivity index can be determined by the future 
observations. As shown by the shaded regions in Figures 3a and 4a, the black hole spin could be constrained in the range of 
0.7 < a. < 0.998 with the emissivity indexes in the ranges of 3.95 < /? < 4.12 and 3.95 < (3 < 4.24 for XTE J1650-564 and 
GRO J1655-40, respectively. 

As pointed out in RM06, HFQPOs are likely to offer the most reliable measurement of BH spins once the correct model 
is known. The observations involved HFQPOs are very complicated, behaving multiple features in different sources, and a 
variety of mechanisms and models are invoked to fit the HFQPOs in BHXBs. The model of the MC process with a nonzero 
torque exerted at the inner edge of accretion disc provides an interpretation for the association of the 3:2 HFQPO pairs with 
the broad Fe K line, and it could be a possible approach to the BH spins of the binaries in the Galaxy. 



4 DISCUSSION 

In this paper, association of the 3:2 HFQPO pairs with the broad Fe K line in XTE J1550-564 and GRO J1655-40 is discussed 
based on MC-II model. It is shown that the 3:2 HFQPO pairs observed in the two sources can be fitted with the observational 
constraints by tuning several model parameters. Some related issues are given as follows. 

(1) It is noted that the mechanisms of creating the two hotspots are different: the inner one is produced predominantly 
by the magnetic torque exerted at the inner edge of the disc, while the outer one arises from the screw instability of the 
large-scale magnetic field. It is helpful to imagine the magnetic field line as an elastic string. The rotating BH twists the field 
line, while the field line tries to untwist itself. Once the toroidal component of the magnetic field is strong enough to satisfy 
the criterion, the screw instability will occurs, just as a twisted elastic string releases its energy under appropriate conditions. 

(2) Some authors argued that the coronal heating in some stars including the Sun is probably related to dissipation of 
currents, and very strong X-ray emissions arise from variation of magnetic fields (Galsgaard & Parnell 2004; Peter et al. 2004). 
Analogously, if the corona exists above the disc in MC-II model, we expect that the corona above Tout might be heated by the 
induced current due to the screw instability of the non-axisymmetric magnetic field. Thus the ratio of the disc bolometric flux 
to the power-law flux is greater for the inner hotspot than that for the outer hotspot, and the feature of X-ray radiation given 
by R02 can be interpreted based on MC-II model with corona, i.e., there exists a systematic shift away from disc dominated 
flux with the increasing power-law flux as the HFQPO pairs shift from the higher frequency to the lower frequency based on 
the observations of XTE J1550-564 and GRO J1655-40. 

(3) Since the two hotspots are produced by two different mechanisms, we can explain the observations that HFQPOs 
do not always appear in simultaneous pairs, because the condition for the magnetic torque and that for the screw instability 
would not be satisfied at the same time. Since the inner and outer hotspots are related to two different mechanisms based 
on the non-axisymmetric magnetic fields in the plunging region and the BH horizon, either the fiuctuations of the magnetic 
fields at different position or the criterion of the screw instability could affect the HFQPO pairs. Unfortunately, we cannot 
discuss this issue in a quantitative way due to the unknown origin of the magnetic field near the BH. 

(4) According to equation (21) the HFQPO frequency is inversely proportional to the BH mass, which providers a strong 
limit to the fits. We can neither fit the HFQPO with 92 Hz in XTE J1550-564 nor the 3:2 HFQPO pair (67, 41Hz) in GRS 
1915-1-105 based on equation (21) for the measured BH masses. Furthermore, the 3:2 HFQPO pair (168, 113Hz) of GRS 
1915-1-105 are not fitted, because a systematic shift away from disc dominated flux with the increasing power-law flux as the 
HFQPO pair shift from the higher frequency to the lower frequency has not been observed in this source. 

(5) As argued above, the higher and lower HFQPO frequencies are related to the inner and outer hotspots rotating 
with the angular velocity located at and rout, respectively. According to equation (21) the HFQPO frequencies are 
inversely proportional to the BH mass, and the spin is related to the HFQPO frequencies via and rout, where the HFQPO 
frequencies are calculated by the angular velocity. Since the evolution time scale of the BH mass and spin are much longer 
than the observation time scale, the two quantities are considered flxed in calculations. 



3:2 HFQPO Pairs and Broad Fe K Line 7 



(6) As pointed out in R02, the 3:2 HFQPO pairs could have a smah, but significant, change in frequency. For example, 
the HFQPO frequencies observed in XTE J1550-564 deviate from 274 Hz and 184 Hz during different outbursts. These small 
deviations in HFQPO frequencies can be calculated by resolving FIE with the adjustable parameters n and S as shown in 
Table 2. 

(7) As argued above, both the non-axisymmetric magnetic field configuration and the corona above the disc surface are 
required for interpreting the 3:2 HFQPO pairs observed in the SPL state of the BHXBs. Thus we can interpret the observation 
that sometimes the 3:2 HFQPO pairs cannot be observed in the SPL state, if the magnetic field configuration is axisymmetric. 
Based on the same reason we predict that the 3:2 HFQPO pairs could be observed with the very steep emissivity index, if the 
magnetic field configuration is non-axisymmetric. However, the HFQPO pairs cannot be produced with the emissivity index 
for the axisymmetric magnetic field configuration. 

(8) According to Bd oc r~" the ratio of the magnetic field at rout to that at ri„ can be written as Rb = {BD)out/{BD)in = 
(rin/rout)" ■ Based on FIE we have the curves of Rb for XTE J1650-564 vary with a, as shown in Figure 5, from which a 
strongly constraint to the BH spin can be found. The BH spin should be very high, or at least greater than some intermediate 
value to avoid the magnetic field at Tout too low to produce an outer hotspot. 

(9) There is less certain evidence that the same 3:2 ratio occurs for QPOs observed in low-mass active galactic nuclei, Sgr 
A*, (Torok 2005a, b; Ashenbach 2004a, b) and in a few nearby Seyferts (Lachowicz et al. 2006). MC-H model can be used to 
fit QPOs pairs observed in the super-massive BH systems, and we obtain the curves of the emissivity index /3, the power-law 
index n and the parameter 5 varying with the BH spin of for Sgr A* as shown in Figure 6. 

We can find from Figure 6 that the parameter n varies from 4.272 to 7.108 with a* decreasing from 0.998 to 0.853 for 
rriBH = 4.4 x lO*", while it varies from 5.040 to 10.239 with a, decreasing from 0.998 to 0.555 for uibh = 2.6 x 10®. 

(10) The parameter n is used to indicate the variation of the large-scale magnetic field with the disc radius, which is a 
key parameter for fitting the 3:2 frequency ratio. Based on FIE we have the variation of the 3:2 frequency ratio with the 
parameter n for the given values of the BH spin as shown in Figure 7. 

Inspecting Figure 7, we find that the frequency ratio decreases monotonically with the increasing param- 
eter n. The greater n corresponds to the less ratio for the given spin, and the greater spin corresponds to 
both the less value and less variation of n for the given ratio. For XTE J1550-564, e.g., the ratio increases 
about 5.7% as the value of n decreases about 10% for a, — 0.5, while the same variation occurs as n only 
decreases about 1.9% for a, = 0.998 . This result imphes that the observed 3:2 frequency ratio is reached for 
a very narrow range of n with the greater BH spin. 

(11) Finally, we discuss the values of the parameters an and am, which are used to indicate the disc thickness at ri„ and 
the strength of the accretion rate, respectively. The emissivity index of XTE J 1550-564 varies with the BH spin for different 
values of a£f and a™ are shown in Figure 8. 

Inspecting Figure 8a, we find that the values of the emissivity index are insensitive to the values of an and am , provided 
that they are limited to an ^0.1 and a™ ^ 1.0. These value ranges are reasonable for thin accretion disc in MC-II model. In 
addition, based on the features of MC-II model, the 3:2 HFQPO pairs are also insensitive to the parameters with the above 
limited values. 
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Figure 3. The curves of the emissivity index /?, the power-law index n and the parameter 6 versus the black hole spin a* for XTE 
J1550-564 are shown in Figures 3a, 3b and 3c, respectively. The solid and dotted lines correspond to the upper and lower limits to the 
black hole mass, respectively. The shaded region in Fig. 3a indicates 3.95 < (3 < 4.12. 



dotted: m_BH=6.0, solid: m_BH=6.6 dotted: m_BH=6.0, solid: m_BH=6.6 dotted: m_BH=6.0, solid: m_BH=6.6 




a* a* 



Figure 4. The corresponding curves for GRO J1655-40 with the same caption as given in Fig.3, except that the shaded region indicates 
3.95 < 13 < 4.24. 



3:2 HFQPO Pairs and Broad Fe K Lme 

Table 1. Calculation results for the 3:2 HFQPO pairs and the emissivity index based on MC-II model. 
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Source 



Input Quantities 



"IBH 



Output Quantities 



XTE J1550-564 276Hz 



184Hz 



0.5 



0.7 



8.4 
10.8 



8.4 
10.8 



0.998 8.4 
10.8 



12.29 
12.60 



2.22 X 10'"' 
4.42 X 10-^ 



10.05 
9.99 



6.37 X 10"'' 
3.29 X 10"* 



6.23 
5.65 



6.22 X 10^^ 
5.52 X 10-3 



4.25 
4.54 



4.12 
4.33 



3.83 
3.95 



GRO J1655-40 



450Hz 



300Hz 



0.5 



0.7 



6.0 
6.6 



6.0 
6.6 



0.998 6.0 



12.44 
12.59 



1.03 X 10"* 
4.60 X 10-5 



9.99 
9.99 



4.45 X lO"" 
3.33 X 10"'' 



5.87 
5.66 



5.79 X 10-3 
5.53 X 10-3 



4.41 
4.53 



4.24 
4.33 



3.90 
3.95 



Notes: The input quantities of XTE J1550-564 and GRO J1655-40 are taken from RM06. The parameter S is defined in equation I I20I I 
as the fraction of the power of the magnetic torque to keep the inner edge at ri„. 



Table 2. Fits of the small changes in HFQPO frequency of XTE J1650-564 with mBH = 10.8. 



Source 



Slower 

(Hz) 
Ratio 



1998-99A 



1998-99B 



2000 Apr 30-May 9 



281.7 ±1.5 
187.8 



277.7 
185.1 ± 3.5 



1.5 



+0.01 
-0.01 



1.5 



-0.03 



269.4 ±2.7 
178.6 

1 t; + 0.02 

-L.tJ I\ Ml 



0.5 



0.7 



0.998 



3-46+o:24 X 4.12 X 10-^ 



12.64 
4.565 



-0.14 
+0.14 

-0.008 
-0.008 



12.60 



+0.57 



4.544 



5.68 



'+D.56 



X 10" 



12.40 



-0.25 
+0.27 



4.502 



+0.014 
-0.013 



-0.06 



9.99 



-0.10 
+0.10 



4 356+"-°°^ 



3.22 X 10-"' 3.57~g ll x 10" 
9.887''-^* 
4.310^ 



9 qs+O-39 
»-»o_0.35 



4.341 



°+0.19 
fO.OlO 



-0.009 



5.47;[;:^i X 10- 



5.51 X 10 



-3 5.59;0:«3 X 10" 



5.61 



-0.02 
+0.03 



3 959+0.003 



5 64+0-05 
3-"*-0.04 



3.951 



n Qqc- + 0.005 

o.yoo_g QQg 



Notes: (1) The super- and sub-scripts in 5, n and /3 correspond to up and down deviations of HFQPO frequencies, respectively. (2) The 
value of & only depends on a* and ^higher i being independent of the deviation of i^jo^er ™ the case of 1998-99B. 
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Figure 6. The curves of the emissivity index f}, the power-law index n and the parameter 5 versus the black hole spin a* for Sgr A* 
are shown in Figures 6a, 6b and 6c, respectively. The solid and dotted lines correspond to the higher and lower limits to the black hole 
mass {ruBH = 4.4 x 10^ and tubh = 2.6 x 10® ). 
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m=8.4, 276Hz fixed 



m=6.0, 450Hz fixed 



1.6 



1 . 55 



o 1-5 



% 1.45 



1.4- 



1.6 



1.55 



o 1.5 - 



10 



10 



Figure 7. The variation of the frequency ratio with the parameter n for (a) XTE J1550-564 and for (b) GRO J1655-40. The solid, 
dashed and dotted Unes correspond to o« = 0.998, 0.7 and 0.5, respectively. 




Figure 8. The curves of einissivity index of XTE J1550-564 versus the black hole spin for (a) am = 0.1, and au = 0.01, 0.1, 0.5 and 
1.0 in solid, dashed, dot-dashed and dotted lines, respectively, and for (b) an = 0.1, and am = 0.01, 0.1, 1.0 and 10.0 in solid, dashed, 
dot-dashed and dotted lines, respectively. The black hole mass is taJsen as tubh = 10.8. 



